Random walk with shrinking steps 
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P. L. Krapivskyi'0 and S. Redneri-Q 

^ Center for BioDynamics, Center for Polymer Studies, 
and Department of Physics, Boston University, Boston, Massachusetts 02215 

We outline the properties of a symmetric random walk in one dimension in which the length of the 
nth step equals A", with A < 1. As the number of steps N ^ oo, the probability that the endpoint 
is at X approaches a limiting distribution P\{x) with many beautiful features. For A < 1/2, the 
support of P\{x) is a Cantor set. For 1/2 < A < 1, there is a countably infinite set of A values for 
which P\{x) is singular, while P\{x) is smooth for almost all other A values. In the most interesting 
case of A = 3 = (V5 - l)/2, Pg{x) is riddled with singularities and is strikingly self-similar. The 
self-similarity is exploited to derive a simple form for the probability measure M{a, b) = Pg{x) dx. 

PACS numbers: 



I. INTRODUCTION 

This article discusses the properties of random walks in 
one dimension in which the length of the n step changes 
systematically with n. That is, at the nth step, a particle 
hops either to the right or to the left with equal proba- 
bility by a distance ,f{n), where ,f{n) is a monotonic but 
otherwise arbitrary function. The usual nearest-neighbor 
random walk is recovered when /(n) — 1 for all n. 

Why should we care about random walks with vari- 
able step lengths? There are several reasons. First, 
for the geometric random walk, where f{n) = A" with 
A < 1, a variety of beautiful and unanticipated features 
arise, 0, 12, 0, 13 illustrated in Fig.^ A very surprising 
feature is that the character of the probability distribu- 
tion of the walk changes dramatically as A is changed 
by a small amount. Most of our discussion will focus 
on this intriguing aspect. We will emphasize the case 
X = g = (a/S — l)/2, the inverse of the golden ratio, 
where the probability distribution has a beautiful self- 
similar appearance. We will show how many important 
features of this distribution can be obtained by exploit- 
ing the self-similarity, as well as the unique numerical 
properties of g. 

There also are a variety of unexpected applications 
of random walks with variable step lengths. One ex- 
ample is spectral line broadening in single-molecule 
spectroscopy. 'sj Here the energy of a chromophore in 
a disordered solid reflects the interactions between the 
chromophore and the molecules of the host solid. For 
a dipolar interaction potential and for two-state host 
molecules with intrinsic energies ±e, the energy of the 
chromophore is proportional to ±er~'^, where rj is 
the separation between the chromophore and the jth 
molecule. This energy is equivalent to the displacement 
of a geometric random walk with f{n) = n~^. 

Another example is the motion of a Brownian particle 
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in a fluid with a linear shear flow, that is, the velocity 
field is Vx{y) = Ay. As a Brownian particle makes pro- 
gressively longer excursions in the zty directions (propor- 
tional to t^/"^), the particle experiences correspondingly 
larger velocities in the x-direction (also proportional to 
t^/^). This gives rise to an effective random walk process 
in the longitudinal direction in which the mean length of 
the nth step grows linearly with n, that is, /(n) = n.0 

Historically, the geometric random walk has been dis- 
cussed, mostly in the mathematics literature, starting 
in the 193O's.0, 0] Interest in such random walks has 
recently revived because of connections with dynamical 
systems. 0, |^ Recent reviews of the geometric random 
walk can be found in Ref. 0; see also Ref. ^3 for a review 
of more general iterated random maps. In contrast, there 
appears to be no mention of the geometric random walk 
in the physics literature, aside from Ref. [Tll |. 

Finally, the geometric random walk provides an in- 
structive set of examples that can be analyzed by classi- 
cal probability theory and statistical physics tools. [l^lT^ 
These examples can serve as a useful pedagogical supple- 
ment to a student's introduction to the theory of random 
walks. As we shall discuss, there are specific values of A 
for which the probability distribution can be calculated 
by elementary methods, while for other A values, there is 
meager progress toward an exact solution, even though 
many tantalizing clues exist. 



II. PICTURE GALLERY 

The displacement of a one-dimensional random walk 
after iV-|- 1 steps has the form xi^ = X]!^=o ^nf{n), where 
each e„ takes on the values ±1 with equal probability. 
Consequently, the mean-square displacement {x%) can be 
expressed as J2n fi''^)'^- this sum is finite, then there is 
a finite mean-square displacement as N —^ 00. In such a 
case, the endpoint probability distribution therefore ap- 
proaches a fixed limit. 14] Henceforth, we will focus on 
the case of geometrically shrinking step lengths, that is, 
/(n) = A" with A < 1. We denote the endpoint proba- 
bility distribution after TV steps by P\{x,N) and limit- 



2 



ing form liniAr^oo P\{x, N) by P\{x). We will show that 
P\{x) exhibits rich behavior as A is varied. 
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FIG. 1: Simulated probability distribution P\{x) for A = 
0.53,0.58,0.61,0.64,0.67, and 0.74 ((a) - (f)). The data for 
each A is based on 10* realizations of 40 steps at spatial res- 
olution 10~^. 



To obtain a qualitative impression of P\{x), a small 
picture gallery of P\{x) for representative values of A > 
1/2 is given in Fig.^ As we shall discuss, when A ~ 1/2, 
^72(2;) — 1/4 for |a;| < 2 and Pi/2(x) = otherwise. 
For A < 1/2, the support of the distribution is fractal 
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(see below). As A increases from 1/2 to approximately 
0.61, P\(x) develops a spiky appearance that changes 
qualitatively from multiple local maxima to multiple lo- 
cal minima (see Fig.^. In spite of this singular appear- 
ance, it has been proven by Solomyakjl^j that the cumu- 
lative distribution is absolutely continuous for almost all 
A > 1/2. On the other hand, Erdos'^ showed that there 
is a countably infinite set of A values, given by the recip- 
rocal of the Pisot numbers in the range (1, 2), for which 
Px{x) is singular. A Pisot number is an algebraic num- 
ber (a root of a polynomial with integer coefficients) all 
of whose conjugates (the other roots of the polynomial) 
are less than one in modulus (more details about Pisot 
numbers can be found in Ref.[la|). It is still unknown, 
however, if these Pisot numbers constitute all of the pos- 
sible A values for which P\{x) is singular. For A > 0.61, 
P\{x) rapidly smooths out and beyond A ^ 0.7, there is 
little visual evidence of spikincss in the distribution at 
the 10~^ resolution scale of Fig. ^ 

There is a simple subset of A values for which singu- 
lar behavior can be expected on intuitive grounds. In 
particular, consider the situations where A satisfies 



JV 



1 - ^ A" = 0. 



This condition can be viewed geometrically as a walk 
whose first step (of length 1) is to the right and N sub- 
sequent steps are to the left such that the walker returns 
exactly to the origin after 7V-|- 1 steps. This positional de- 
generacy, in which points are reached by different walks 
with the same number of steps, appears to underlie the 
singularities in P\{x). The roots of Eq. give the solu- 
tion A = 5 EE (^/5 - l)/2 w 0.618, 0.5437, 0.5188, 0.5087, 
0.5041 ... for = 2, 3, ... As we shall discuss, the largest 
in this sequence, X = g (the inverse of the golden ratio), 
leads to especially appealing behavior and the distribu- 
tion Pg{x) has a beautiful self-simila rity as well as an 
infinite set of singularities. |1 HItL HtL ITsj 



III. GENERAL FEATURES OF THE 
PROBABILITY DISTRIBUTION 

For A < 1/2, the support of P\{x), namely the subset 
of the real line where the distribution is non-zero, is a 
Cantor set. To understand how this phenomenon arises, 
suppose that the first step is to the right. The maximum 
displacement of the subsequent walk is A/(l — A). Conse- 
quently, the endpoint of the random walk necessarily lies 
within the region [l — j^j, 1 + y^j] ■ Thus the support of 
P\{x) divides into two non-overlapping regions after one 
step (see Fig.|21). This same type of bifurcation occurs at 
each step, but at a progressively finer distance scale so 
that the support of P\{x) breaks up into a Cantor set. 

The evolution of the support of the probability distri- 
bution can be determined in a precise way by recasting 





FIG. 2: Fragmentation of the support of the probability dis- 
tribution during the first two steps of the shrinking random 
walk when A < 1/2. After k steps, the support fragments into 
2^ intervals each of length 2A'=/(1 - A). 



the geometric random walk into the random map, 

x' = ±l + \x, (1) 

which describes how the position of a particle changes in 
a single step. We can see that this map is equivalent to 
the original random walk process by substituting x' for x 
on the right-hand side of Eq. Q and iterating. Therefore 
the map is equivalent to the random sum x — ^nX^ ■ 
Similarly, the equation for the probability distribution 
at the A^th step satisfies the recursion 

relationjUm 

P^[x,N) = \\p^{^,N-l)+P>^{^^,N-l) 



A 



(2) 

Equation jSJ states that for a particle to be at x at the 
A^th step, the particle must either have been either at 
(a;-|-l)/A or at {x—l)/\ at the previous step; the particle 
then hops to x with probability 1/2. For any A < 1, the 
probability distribution necessarily approaches a fixed 
limit P\{x) for N — > oo. Consequently, P\{x) remains 
invariant under the mapping given by Eq. J^l and thus 
satisfies the invariance condition 



1. 



A 



Pa 



Is 



A 



(3) 



Because P\{x) can have a singular appearance, it often 
is more useful to characterize this distribution by the 
probability measure Mx{a,b) defined by 



(4) 



Mxia,b) = / da; Pa (a;). 



The integral in Eq. Q smooths out singularities in Pa 
itself, and we shall see that it is more amenable to theo- 
retical analysis. The invariance condition of Eq. Q can 
be rewritten in terms of this measure as 



2AfA(a,6) = A/a 



15-1 



Ah 



1 b+1^ 



.(5) 
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Equation (jSJ can now be used to determine the support 
of M\. Clearly, the support of M\ lies within the interval 
J\ = [-a;max,a;max], with Xmax = 1/(1 -A). For A < 1/2, 
the map ^ transforms the interval J\ into the union of 
the two non-overlapping subintervals (see Fig. 12), 



(l-2A)i r(l-2A) 



(1-A)' (l-A)J' L(l-A)'(l-A) 



(6) 



By restricting the map ^ to these two subintervals, we 
find that they are transformed into four non-overlapping 
subintervals after another iteration. If we continue these 
iterations ad infinitum, we obtain a support for M\ that 
consists of a collection of disjoint sets that ultimately 
comprise a Cantor set.0 

On the other hand, for A > 1/2, the map again trans- 
forms J\ into the two subintervals given in Eq. ©, but 
now these subintervals are overlapping. Thus the support 
of Pa fills the entire range [-Xmax, a;inax]- 



IV. EXACT DISTRIBUTION FOR A = 2"^''™ 

In this section, we derive Pa by Fourier transform 
methods for A = 2"!/™. We wiU illustrate that the dif- 
ferent values of m turn out to be exactly soluble because 
of a set of fortuitous cancellations in the product form 
for the Fourier transform of the probability distribution. 

For a general random walk process, the probabil- 
ity P{x, N) that the endpoint of the walk is at x at 
the {N + 1 ) st ep obeys the fundamental convolution 
eauation|ia [l3 | 



P(x, N)^Y. -x',N- 1)pn{x'), 



(7) 



where Pn{x) is the probability of hopping a distance x 
at the N*'^ step. Equation {Tj) expresses the fact that to 
reach x after {N + 1) steps the walk must first reach a 
neighboring point x~x' after iV steps and then hop from 
x — x' to X at step (N -f 1). The convolution structure of 
Eq. Q cries out for employing Fourier transforms. Thus 
we introduce 



PN{k) 
P{k,N) 



Pn{x) e'''^ dx 

o 

P{x, N) e''^^ dx, (8) 



and substitute these forms into Eq. ((T)). If the random 
walk is restricted to integer-valued lattice points, these 
integrals become discrete sums. The Fourier transform 
turns the convolution in x into a product in fc-space.|l9| 
and therefore Eq. ||7J| becomes the recursion relation 

P(k,N) = P{k,N -l)pN{k). (9) 

We now iterate Eq. ^ to obtain the formal solution 



N 



P{k,N)^ P{k,0) Y[pn{k). 



(10) 



Generally, we consider the situation where the random 
walk begins at the origin. Thus P(a;,0) = S^fi and cor- 
respondingly from the second of Eq. (|SJ), P(fc, 0) — 1. To 
calculate P{x,N), we evaluate the product in Eq. 110|) 
and then invert the Fourier transform. (IZJ To simplify 
the notation for the examples of this section, we define 
^m{x) — P2-i/m{x). We explicitly consider the cases 
TO = 1, 2, and 3. From these results, the qualitative 
behavior for general to is then easy to understand. 
The single-step probability distribution at the nth step 

is 

p„(x) = l-[5{x-\") + Six + X")], 



and the corresponding Fourier transform is 

Pn{k) = / p„(a;)e*''^ = cos(fcA"). 



(11) 



The Fourier transform of the probability distribution af- 
ter N+1 steps is the product of the Fourier transforms of 
the single-step probabilities. 0,^3 Thus from Eqs. H10|l 
and Hll|l we have 



N 



Px{k,N) = Yl cos(/cA"). 



(12) 



n=0 



We now apply this exact solution to the illustrative 
cases of A = 2^^"^. For the simplest case of A = 2~^, the 
step length systematically decreases by a factor of 2. By 
enumerating all walks of a small number of steps, it is 
easy to see that the probability distribution is uniformly 
distributed on a periodic lattice whose spacing shrinks 
by a factor of two at each step. This is a precursor of 
the uniform nature of ni(a;) in the N —^ oo limit. Al- 
gebraically, the product in Eq. (|12|) can be simplified by 
using the trigonometric half-angle formula to yield 



Iii{k,N) 



cos k cos(fc/2) . . 
sin(2A:) sinfc 



2sinfc 2sin(fc/2) 

sin(2fc) 
2^+isin(A:/2^) 
sin(2fc) 
2k 



cos(fc/2^) 

sin(fc/2^-i) 



2sin(fc/2^) 



as — > oo. 



(13) 



Thus the inverse Fourier transform gives an amazingly 
simple result. The probability distribution is merely a 
square- wave pulse, with Hi = 1/4 on the interval [—2, 2] 
and Hi = otherwise. 

The distribution for A — 2^^!"^ can be calculated sim- 
ilarly. The successive cancellation of adjacent numera- 
tors and denominators as in Eq. (|13|1 still occurs. These 
cancellations become more apparent by separating the 
factors that involve sin(A;/2^) and sin(A:/2-'"'"^/^). Then 
by following exactly the same steps as those leading to 
Eq. H13|l . we obtain the Fourier transform 



n=0 



n2(fc) 



sin(2fc) sin(V2fc) 



(14) 
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This product form has a simple interpretation in real 
space. If we partition the walk into odd steps (1,3,5,...) 
and even steps (2,4,6, . . .), then both the odd and even 
steps are separately geometrical random walks with A = 
2~^, but with the initial step length of the odd walk 
equal to 1 and that of the even walk equal to l/\/2- In 
real space, the full probability distribution for A = 2^^/^ 
is just the convolutions of the constituent distributions 
for these odd and even walks. Thus in Fourier space, the 
full distribution is just the product of the constituent 
distributions, as given in Eq. H14|) . 



Oh 




FIG. 3; Probability distributions for the GW for A = 2-1/™ 
for m — 1, 2, and 3 (dotted, dashed, and solid, respectively). 



To invert the Fourier transform in Eq. H14(l by a direct 
approach is straightforward but unwieldy, and the details 
are given in Appendix A. Another approach is to use the 
fact that the probability distribution is the convolution of 
two rectangular-shaped distributions - one in the range 
[—2,2] for the odd-step walks and the other in [— v^, V^] 
for the even-step walks. Thus 



n2(x) 



: f dx'Uiix') X V2Ui(- 

J — OQ ^ 



V2 



(15) 



Either by this direct approach or by straightforward 
Fourier inversion as given in the Appendix, the final re- 
sult is (Fig. EJ) 



Tl2{x) 



2+V2> 



{\x\ <2-V2) 

(2-x/2< \x\ <2 + V2)(16) 
(|a;| >2 + ^/i) 



Thus the distribution is continuous, but its first deriva- 
tive is discontinuous at the four points x ~ ±2 ± \/2. 

Continuing this same train of logic, the solution for 
general A = 2-^/™ is 



n„(fc) - 



njLi sin(2^/'"fc) 
2(™+i)/2fc™ ' 



(17) 



For example, for A = 2 '^/'^, the resulting probability 
distribution in real space is 



64 (^4 ^\ 



= < 



64 (^4 ■ 
64(^4 ■ 



x"^) {\x\<Xi) 

2xxi) (xi < \x\ < X2) 



xl - 2x{xi + X2)) 



4^{x — X4y 



where 



Xl 

X4 



{X2 < \x\ < X3) 
{X3 < \x\ < X4:) 
{\x\ > X4) 



-2 + 2^/^ + 2^/^ 
+2 - 2^/^ + 2^/^ 
+2 + 2^/^ - 2^/^ 
+2 + 2^/^ + 2^/^ 



0.08473 
1.6725 
2.3275 
4.8473 



This distribution contains both linear and quadratic seg- 
ments such that the first two derivatives of II3 are con- 
tinuous, but the second derivative is discontinuous at the 
joining points Xj, for j = 1, 2, 3, 4 (Fig.lSJ. Generally, for 
A — 2"-^/™, the distribution has continuous derivatives 
up to order to — 1 , while the TOth derivative is discontin- 
uous at 2m points. As to ^ cx3, the distribution progres- 
sively becomes smoother and ultimately approaches the 
Gaussian of the nearest-neighbor random walk. 

A final note about the Fourier transform method is 
that it provides a convenient way to calculate the mo- 
ments, 



{x'^) = / x'''Px{x)dx, 



(19) 



for all values of A.|20| By expanding Px{k) is a power 
series for small fc, we have 



Px{k) 



dx P\ [x) e 



ikx 



dx Px{x) (1 + ikx - k^x^/2\ + . . .) 



= 1-^ + ^ + .... (20) 

That is, the Fourier transform contains all the moments 
of the distribution. For this reason, the Fourier transform 
is often termed the moment generating function. 

We take the Fourier transform of the probability dis- 
tribution of the geometric random walk and expand this 
expression in a power series in k to give 

Px{k) — cosfc cos(Afc) cos(A^fc) . . . 



1 



fc2 

' 2! 
fc2 



1- 



2! 



[1- 



(A2fc)2 



2! 



l-y(l + A2 + A4 + ...) + 0(fc4). 



If we equate the two power series (|20f) and (|21f) term by 
term, we obtain 



1 



{x') 



1 - A2 

1 

1~A4 



6A' \ 
1~A2J' 



(21a) 
(21b) 
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Moments of any order can be obtained by this approach. 



V. GOLDEN WALK 

Particularly beautiful behavior of Pxix) occurs when 
\ — g (see Fig. ^d)). Unfortunately, a straightforward 
simulation of the geometric random walk is not a practi- 
cal way to visualize fine-scale details of this probability 
distribution accurately because the resolution is neces- 
sarily limited by the width of the bin used to store the 
distribution. We now describe an enumeration approach 
that is exact up to the number of steps in the walk. 



A. Enumeration 



Xiq to a first-order polynomial. To this end, we succes- 
sively use the defining equation = 1 — g to reduce 
all powers of g to first order. When we apply this re- 
duction to g", we obtain the remarkably simple formula 
g^ = (— — gFn), where F„ is the nth Fibonacci 
number (defined by Fn — Fn-i + -Fn-2 for n > 2, with 
Fi = F2 = 1). For the golden walk, we now use this con- 
struction to reduce the location of each endpoint, which 
is of the form CnS", to the much simpler form A+Bg, 
where A and B are integers whose values depend on the 
walk. By this approach, each endpoint location is ob- 
tained with perfect accuracy. The resulting distribution, 
based on enumerating the exact probability distribution 
for N < 29, is shown in Fig. [S] at various spatial resolu- 
tions. At iV = 29 this distribution is exact to a resolution 
of 10-^. 



It is simple to specify all walks of a given number of 
steps. Each walk can be represented as a string of bi- 
nary digits with representing a step to the left and 1 
representing a step to the right. Thus we merely need 
to list all possible binary numbers with a given number 
of digits N to enumerate all A^-step walks. However, we 
need a method that provides the endpoint location with- 
out any loss of accuracy to resolve the fine details of the 
probability distribution. 




FIG. 4: First 5 steps of the golden walk enumeration tree. 
Notice that the distance between adjacent endpoints can only 
be either 2g" or 2g"-\ 

The basic problem of determining the endpoint loca- 
tions may be illustrated by enumerating the first few 
steps of the walk and representing them as the branches 
of a tree as shown in Fig. 01 For X = g, neighboring 
branches of the tree may rejoin because of the existence 
of 3-step walks (one step right followed by two steps left 
or vice versa) that form closed loops. For large N, the 
accuracy in the position of a walk is necessarily lost by 
roundoff errors if we attempt to evaluate the sum for the 
endpoint location, = J2n ^^g^, directly. Thus the 
recombination of branches in the tree will eventually be 
missed, leading to fine-scale inaccuracy in the probability 
distribution. 

However, we may take advantage of the algebra of 
the golden ratio to reduce the iVth-order polynomial in 



B. Self-Similarity 

Perhaps the most striking feature of the endpoint dis- 
tribution is its self-similarity, as sketched in Fig. |B| No- 
tice, for example, that the portion of the distribution 
within the zeroth subinterval J'^ = [— g] is a micro- 
cosm of the complete distribution in the entire interval 
J = [~g~^ , g~^]- In fact, we shall see that the distri- 
bution within J° reproduces the full distribution after 
rescaling the length by a factor g~^ and the probability 
by a factor of 3. Similarly, the distribution in the first 
subinterval = [<7j 1 + <7^] reproduces the full distribu- 
tion after translation to the left by 1, rescaling the length 
by <?~^, and rescaling the probability by 6. A similar con- 
struction applies for general subintervals. 

To develop this self-similarity, it is instructive to con- 
struct the symmetries of the probability distribution. 
Obviously, Pg{x) is an even function of x. That is, 

Pg{x) = P,{-~x). (22) 

In fact, there is an infinite sequence of higher-order sym- 
metries that arise from the evenness of Pg{x) about the 
endpoints after 1 step, 2 steps, 3 steps, etc. 

For example, the first higher-order symmetry is 

Pgil + X) ^ Pg{l - X) (23) 

for \x\ < g^. Equation (|23|l expresses the symmetry of 
Pg {x) about x = 1 for the subset of walks whose first step 
is to the right. We can ignore walks whose first step is 
to the left because the rightmost position of such walks 
is — 1 + 5 + .g^ + . . . = g = 1 — 5^- Thus within a distance 
of g^ from x ~ 1, only walks with the first step to the 
right contribute to the distribution within this restricted 
range. The probability distribution must therefore be 
symmetric about x = 1 within this same range. 
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FIG. 5: Probability distribution of the golden walk for a 29- 
step enumeration at spatial resolution 10~^, 10""^ and 10"* 
(a) - (c), respectively. In (c), the line joining successive points 
is not shown so that details of the distribution remain visible. 



Continuing this construction, there are infinitely many 



symmetries of the form 



^9(E5" + ^)-^9(E-9"-^)' (24) 

n=0 n=0 

with k = 1,2,... that represent reflection symmetry 
about the point that is reached when the first k steps are 
all in one direction. The kth symmetry applies within 
the range \x\ < g''^^ ■ 

We now exploit these symmetries to obtain a simple 
picture for the measure of the probability distribution, 
Mg. We start by decomposing the full support J into 
the contiguous subintervals that span the successive lobes 
of the distribution, as shown in Fig. We label these 
subintervals as J° = (— g, g), = (1 — g^, 1 + g^), = 



(l + 5-5',l 



g ), etc.; there are also mirror image 



intervals to the left of the origin, J^'^ = —J^ . 

We now use the invariance condition of Eq. Q to de- 
termine the measures of these fundamental intervals . 
For = i~g,g), this invariance condition yields 



Mg(~g,g) 



[Mg{-{2 + g),~g) + Mg{g,2 + g)] 

(25) 



= ^[l-Mg{-g,g)], 



where the second line follows because of left-right sym- 
metry and because the intervals {—{2 + g), —g), {g,2 + g), 
and {—g,g) comprise the entire support of a normal- 
ized distribution. We therefore obtain the remarkably 
simple result that the measure of the central interval is 
Mg{~9:g) = 1/3- If we apply the same invariance con- 
dition to J\ we find Mg{J^) = ^Mg{J°). By continuing 
this same reasoning to higher-order intervals , we find, 
in general, that the measure of the fcth interval is one- 
half that of the previous interval (see Fig. El). Thus we 
obtain the striking result 



1 



3-2l'^l 



(26) 



C. Singularities 

Another intriguing feature of Pg{x) is the existence of 
a series of deep minima in the distribution. Consider, 
for example, the most prominent minima at x = ±g (see 
Fig. O . The mechanism for these minima is essentially 
the same reason that g is sometimes termed the "most 
irrational" real number - that is, most difficult to ap- 
proximate by a rational number. ,21j In fact, there is only 
a single trajectory of the random walk in which the end- 
point of the walk reaches namely, the trajectory that 
consists of alternating the steps, 1 — g + g^ — g^ + g'^ — . . .. 
If there is any deviation from this alternating pattern, 
the endpoint of the walk will necessarily be a finite dis- 
tance away from x — g. This dearth of trajectories with 
endpoints close to a; = g is responsible for the sharp min- 
imum in the probability distribution. 
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FIG. 6: Sketch to iUustrate the symmetry and self-similarity 
of Pg{x). The dashed curve is the probability distribution 
when the first step is to the right. The full probability distri- 
bution is the sum of the dashed curve and an identical (but 
shifted) curve that stems from the distribution when the first 
step is to the left. The measures associated with each lobe 
of P\{x) (top) and the spatial extent of each lobe (bottom) 
are indicated. Notice that the left extreme of the restricted 
distribution coincides with the first minimum of the full dis- 
tribution. 

More generally, this same mechanism underlies each of 
the minima in the distribution, including the singularity 
as a: ^ aJinax- For each local minimum in the distribution, 
the first n steps of the walk must be prescribed for the 
endpoint to be within a distance of the order of to the 
singularity. However, specifying the first n steps means 
that the probability for such walks can be no greater 
than 2~". It is this reduction factor that leads to all the 
minima in the distribution. 

For simplicity, we focus on the extreme point in the 
following; the argument for all the singularities is similar. 
If the first n steps are to the right, then the maximum 
distance A between the endpoint of the walk and Xmax 
arises if the remaining steps are all to the left. Therefore 

A™ /111 I 'n\ I n+1 , n+2 i o ^ 

= a;,nax -(1+5 + .. .+5 )+5 +.9 + ■ ■ ■ = ^9 

Correspondingly, the total probability to have a random 
walk whose endpoint is within this range is simply 2~". 

For X near Xmax, we make the fundamental assump- 
tion that P\{x) ~ {xraa^ — xY ■ Although this hypothesis 
appears difficult to justify rigorously for general values 
of A, such a power law behavior arises for A = 2"-'^/™, 
as discussed in Sec. III. We assume that power-law be- 
havior continues to hold for general A values. With this 
assumption, the measure for a random walk to be within 
the range A = Xmax — a; of aJmax is M (A) ^ A^+'^. How- 
ever, because such walks have the first n steps to the 
right, Af(A) also equals 2~". If we write In A/ = — nln2. 



In A = +(n — l)ln5 + ln2, and eliminate n from these 
relations, we obtain M(A) ~ A''^^/ in(i/s) f^^s,\\y^ 

Pg(A) - A"i+'"2/'"(i/f). (27) 

This power law also occurs at each of the singular points 
of the distribution because the underlying mechanism is 
the same as that for the extreme points. 

The same reasoning applies mutatis mutandis near the 
extreme points for general A, leading to the asymptotic 
behavior Pa (A) ^ i^-i+\n2/ _ particular, this 
reason gives, for the tail of H^, the limiting behavior 
A™"^, in agreement with the exact calculation in Sec. HI. 



VI. DISCUSSION 

We have outlined a number of appealing properties 
of random walks with geometrically shrinking steps, in 
which the length of the nth step equals A" with A < 1. A 
very striking feature is that the probability distribution 
of this walk depends sensitively on the shrinkage factor 
A and much effort has been devoted to quantifying the 
probability distribution. We worked out the probability 
distribution for the special cases A — 2~^/™ where a solu- 
tion is possible by elementary means. We also highlighted 
the beautiful self-similarity of the probability distribu- 
tion when A = {sf^ — l)/2. Here, the unique features 
of this number facilitate a numerically exact enumera- 
tion method and also lead to very simple results for the 
probability measure. 

We close with some suggestions for future work. What 
is the effect of a bias on the limiting probability distribu- 
tion of the walk? For example, suppose that steps to the 
left and right occur independently and with probabilities 
p and 1 — p respectively. It has been proven 22] that the 
probability distribution is singular for A < p^(l — p)^^^ 
and is continuous for almost all larger values of A. This 
is the analog of the transition at A = 1/2 for the isotropic 
case. What other mysteries lurk within the anisotropic 
system? 

Are there interesting first-passage characteristics? For 
example, what is the probability that a walk, whose first 
step is to the right, never enters the region x < by the 
1 ?ith step? Such questions are of fundamental importance 
in the classical theory of random walks, 13] and it may 
prove fruitful to extend these considerations to geometric 
walks. Clearly, for A < 1, this survival probability will 
approach a non-zero value as the number of steps TV — > 
oo. How does the survival probability converge to this 
limiting behavior as a function of the number of steps? 
Are there qualitative changes in behavior as A is varied? 

What happens in higher spatial dimensions? This ex- 
tension was suggested to us by M. Bazant.[2^ There are 
two natural alternatives that appear to be unexplored. 
One natural way to construct the geometric random walk 
in higher dimensions is to allow the direction of each step 
to be isotropically distributed, but with the length of the 
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nth step again equal to A". Clearly, if A <C 1, the proba- 
bility distribution is concentrated within a spherical shell 
of radius 1 and thickness of the order of A/ (1 — A). As 
A is increased, the probability distribution eventually de- 
velops a peak near the origin. What is the nature of 
this qualitative change in the probability distribution? 
Another possibility (2^ is to require that the steps are 
always aligned along the coordinate axes. Then for suf- 
ficiently small A the support of the walk would consist 
of a disconnected set, while as A increases beyond 1/2, 
a sequence of transitions similar to those found in one 
dimension may arise. 
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Similarly, the contributions of the remaining three inte- 
grals are —Tr{x2+x), —^{xi — a;),and —Tr{xi +x), respec- 
tively. As a result, we find, for < x < xi, 

{X2- X) + {X2 +X)- {Xl - X) - {Xl +x) 1 

n^(^) ^ ¥n = i- 

Finally, for xi < x < X2, we must close the contour in the 
upper half-plane for the first two terms in Eq. ljA2|l and in 
the lower-half plane for the latter two terms. Evaluating 
each of the residues, we now obtain, for xi < x < X2 

n.w4{i-i-)^ 



APPENDIX A: FOURIER INVERSION OF THE 
PROBABILITY DISTRIBUTION FOR A = T^'^ 

For A = 2^^/^, we write Px{k) in the form 
sin(2fc) sin(V2fc) 



Tl2{k) 



23/2/,2 



2 



2V2/t 

where xi = 2 — \/2 and X2 ^ 2 + \/2. The inverse Fourier 
transform is 

U2{x) = — / n2(fc) e-^'=" dk. (A2) 

To evaluate the integral, we extend it into the complex 
plane by including a semi-circle at infinity. The outcome 
of this inverse transform depends on the relation between 
X and Xl, X2- 

For X > X2, we must close the contour in the lower 
half-plane for each term, so that the semi-circle contri- 
bution is zero. We must also indent the contour around 
an infinitesimal semi-circle about the origin to avoid the 
singularity at fc = 0. If x > X2, the residues associated 
with each term in the integrand cancel, and we obtain 
n2(x) = 0. 

For < X < Xl, we must close the contours in the 
upper half-plane for the first and third terms, and in 
the lower-half plane for the complementary terms. The 
contribution of the first integral is proportional to 



e 



fc2 



■ dk — iTrRes 



ik{x2—3;) 



/c2 



= -7I'(X2 

k=0 
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